Line Polarization of Molecular Lines at Radio Frequencies: The 

case of DR21(OH) 



P. C. Cortes^ 
and 

in '■ R. M. Crutcher^ 

o : 

^ ! and 
9^. W. D. Watson^ 

<: 

(N ■ 

^ : ABSTRACT 
> ■ 

OO I We present polarization observations in DR21(0H) from thermal dust emis- 

^ I sion at 3 mm and from CO J=l line emission. The observations were ob- 

I tained using the Berkeley-IUinois-Maryland Association array. Lai et al. (2003) 

[ observed this region at 1.3 mm for the polarized continuum emission, and also 

^ I measured the CO J = 2 ^ 1 polarization. Our continuum polarization results 

'-^1 are consistent with those of Lai et al. (2003). However, the direction of the linear 

^ I polarization for the J = 1 — > is perpendicular to that of the CO J = 2 ^ 1 po- 

?-j I larization. This unexpected result was explored by obtaining numerical solutions 

^ I to the multilevel, radiative transfer equations for a gas with anisotropic optical 

I depths. We find that in addition to the anisotropic optical depths, anisotropic 

I excitation due to a source of radiation that is external to the CO is needed to 

^ ' understand the orthogonality in the directions of polarization. The continuum 

■ - - ' emission by dust grains at the core of DR21(0H) is sufficient to provide this 

external radiation. The CO polarization must arise in relatively low density ( 
~ 100 cm~^ ) envelope gas. We infer B ~ lO/iG in this gas, which implies 
that the envelope is subcritical. 
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1. Introduction 

The star formation process is one of the most comphcated problems in current astro- 
physical research. Its study includes several physical parameters, of which the magnetic 
field is the least observed. Magnetic field observations are divided into measurements of 
the Zeeman effect (in order to obtain the magnetic field strength in the line of sight), and 
linear polarization observations of star forming regions. Dust polarization is believed to be 
perpendicular to the magnetic field under most conditions (Lazarian 2003); hence, dust 
polarization has been used as a major probe for the magnetic field geometry. In order to 
efficiently measure the dust polarization and infer information about the magnetic field mor- 
phology, high resolution observations are required. The BIMA millimeter interferometer has 
been used to obtain high-resolution maps in several star forming cores (Rao et al. 1998; 
Girart & Crutcher 1999; Lai et al. 2001, 2002, 2003). These results show fairly uniform 
polarization morphologies over the main continuum sources, suggesting that magnetic fields 
are strong, and therefore, should not be ignored in star formation theory. 

The linear polarization of spectral line radiation from the molecular gas has been sug- 
gested to occur under anisotropic conditions (Goldreich & Kylafis 1981a). The prediction 
is that linearly polarized radiation of a few percent is likely to be detected from molecular 
clouds and circumstellar envelopes in the presence of a magnetic field. The direction of this 
polarization should be either parallel or perpendicular to the magnetic field, depending on 
the angles between the line of sight, the magnetic field, and the direction associated with the 
anisotropic excitation (Goldreich & Kylafis 1981b). This process is known as the Goldreich 
- Kylafis effect. Combining this effect with observation of the polarized emission from dust 
grains, we can probe the magnetic field morphology in the plane of the sky. 

We observed the young and massive star forming region DR21(0H). We measured dust 
polarization at 3 mm and the CO J = 1 — line polarization. These observations were 
compared with the Lai et al. (2003) polarization observations at 1 mm and CO J = 2 ^ 1, 
particularly the line polarization that has more extended emission than the dust. We found 
a 90° difference in the position angle when comparing the line polarization data for the two 
transitions. This difference motivated a numerical study to explore polarized line emission 
orientation in different transitions. 

The paper is divided into three major sections. In Section 2 we discuss the source 
and observation procedure, the line and dust polarization observations, and the comparison 
with the Lai et al. (2003) results. The numerical calculational procedure and results are 
described in Section 3. Section 4 contains the discussion and summary. 
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2. Observations and data reduction 

2.1. Source description 

DR21(0H) is a young massive star forming region located in the -3 km DR21/W75S 
molecular cloud complex, centered at RA: 20:39:00.7 and DEC: 42:22:46.7 (J2000 coordi- 
nates). This region is also at the north-eastern part of the giant Cygnus-X H II complex. 

The distance to DR21(0H) is assumed to be 3 kpc; however, the value is uncertain. Dickel 
et al. (1978) used a value of 2 kpc. From thermal emission of dust Woody et al. (1989) 
resolved two compact cores in DR21(0H), MMl and MM2, with a total mass of ~ 125 Mq. 
The MMl component is the brighter one, with an integrated flux of 0.27 Jy at 2.72 mm 
(Mangum et al. 1991). DR21(0H) has been extensively mapped in CO by Dickel et al. 
(1978); Woody et al. (1989); Mangum et al (1991); Chandler et al. (1993), and in CS by 
Plambeck & Menten (1990) and Chandler et al. (1993). DR21(0H) is also known for its as- 
sociation with maser emission from OH (Norris et al. 1982), H2O (Genzel & Downes 1977) 
and CH3OH (Batrla & Menten 1988; Plambeck & Menten 1990). DR21(0H) has also been 
observed in the far-infrared (Williams et al. 1974; Harvey et al. 1986) and sub-millimeter 
(Gear et al. 1988; Richardson et al. 1994). No centimeter-wavelength continuum sources 
have been observed in DR21(0H) (Johnston et al. 1984); therefore, Hll regions have not 
yet developed, so it appears to be in an early stage of evolution. This makes DR21(0H) a 
good candidate for early star formation and magnetic field studies. 

Magnetic field observations have been carried out measuring Zccman splittings. The 
CN line Zeeman splitting has been detected in both cores (the CN line traces n^/j ~ 10^-10^ 
cm~^) giving a line-of-sight magnetic field strength. Bios of -0.4 mG for MMl and -0.7 mG for 
MM2 (Crutcher 1999a). Lai et al. (2003) observed the CO J = 2 ^ 1 molecular line and 
the 1.3 mm dust continuum simultaneously using the BIMA array. They obtained a detailed 
polarization map for the whole region. Their results show a remarkably uniform polarization 
pattern for the line over the main two continuum sources, while the dust polarization appears 
to be mostly perpendicular to the line. These results seem to be consistent with theoretical 
predictions. 

2.2. Observation procedure 

We observed DR21(0H) from October to November 2001, mapping the continuum emis- 
sion at 3 mm and the CO J = 1 ^ molecular line (at 115.2712 GHz), using the BIMA 
array in C configuration. We set the digital correlator in mode 8 to observe both the con- 
tinuum and the CO J = 1 ^ line simultaneously. The 750 MHz lower side band was 
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combined with 700 MHz from the upper side band to map the continuum emission, leaving 
a 50 MHz window for the CO hne observation (at a resolution of 2.06 km s~^). Each BIMA 
telescope has a single receiver and thus the two polarizations must be observed sequentially. 
A quarter wave plate to select either R or L circular polarization is alternately switched in to 
the signal path ahead of the receiver. Switching between polarizations was sufficiently rapid 
(every 11.5 seconds) to give essentially identical uv-coverage. Cross-correlating the right (R) 
and left (L) circularly polarized signals from the sky gave RR, LL, LR, and RL for each 
interferometer baseline, from which maps in the four Stokes parameters were produced. The 
instrumental polarization was calibrated by observing 3c279, and the "leakages" solutions 
were calculated from these observations. The phase calibrator used was BL Lac. We used 
the same cahbration procedure described by Lai et al. (2001). 

The Stokes images I, U, Q and V were obtained by Fourier transforming the visibihty 
data using natural weighting. The 3 mm continuum synthesized beam had a major axis of 
8.4" and a minor axis of 6.9", with a position angle (RA.) of 35°. The MIRIAD (Sault & 
Kileen 1995) package was used for data reduction. 

2.3. 3 mm Continuum results 

Our 3 mm continuum results did not resolve the two continuum sources (MMl and MM2) 
resolved by Lai et al. (2003). The larger beam at 3 mm gives poorer resolution than the 1.3 

mm observations (with BIMA at the same configuration). Our 3 mm continuum observations 
have a peak emission of 0.22 Jy beam'^ centered at 20''39™00^7 RA and 42°22'46".7 DEC 
in J2000 coordinates. This center coincides roughly with the center of the MMl continuum 
source. Johnston et al. (1984) found no continuum emission at 14.5 GHz from DR21(0H). 
Therefore we assume that free-free emission is negligible and all the continuum radiation 
comes from dust emission. The peak emission is consistent with Mangum et al. (1991), who 
obtained 0.192 Jy beam~^ at 2.7 mm with a comparable beam size. 

Polarization detection is very sensitive to the signal to noise ratio. Because of bad 
weather conditions during part of our observations we did not achieve the same level of 
polarization sensitivity as Lai et al. (2003) did. We have 3a polarization detections (figure 
1) scattered over 3 positions in the map. The dust polarization detected is consistent with 
the result of Lai et al. (2003). 
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2.4. CO J —1 ^ Observational results 

Figure 2 shows the Stokes I spectra integrated over a region that contains the MMl 
and MM2 sources (Lai et ah 2003). This region covers the box from RA: 20'*39™2* - 
20^38"59^5 to DEC:42°22'35" - 42°23'2". Plotted over the Stokes I are the P.A. and the 
fractional polarization at each velocity channel. The Stokes I spectra generally agree with 
the CO J = 2 — > 1 spectra from Lai et al. (2003). In our case the spectrum has a peak 
at visr — —8 km s~^, a minimum at visr — —2 km s~^, and a second peak at visr — 2 km 
s'"''". The velocity of these peaks are different from previous CS and CN observations, which 
are at -5 km s~^ (Chandler et al. 1993; Richardson et al. 1994), and -5 km s^^ and -1 
km s^^ (Crutcher 1999a). The CO line is optically thick (Figure 2) and the dip, or zero 
emission range, observed in both CO J = 2 ^ 1 and J = 1 ^ lines, is probably due 
to self-absorption. In CO J = 1 — > 0, the dip covers a velocity range from -6 km s~^ to 2 
km s~^ which coincides with the C^^S J = 2 — > 1 emission peak range from Chandler et al. 
(1993), which goes from -6 km s~^ to 1.5 km s~^. In the same velocity range ^^CO emission 
was detected (Dickel et al. 1978) from -3 km s~^ to 0.5 km s~^. C^^S traces higher densities 
than CO, and ^^CO is more optically thin than CO. The presence of peak emission from 
these lines in the same velocity range of our dip suggests that our CO observations may 
come from a different region from the higher density gas (which probably comes from the 
core). Additionally, it is expected that CO polarized line emission will arise from optically 
thin regions (r ~ 1) (Goldreich & Kylafis 1981a; Deguchi & Watson 1984). Therefore, it 
is likely that the CO polarized emission that we detected comes from an envelope around 
DR21(0H). 

Figure 3 shows our CO polarization map at Visj. =-8 km s^^. This map represents the 
highest significance polarization distribution in our observations. It also coincides with the 
peak emission of the CO. Figure 4 shows our polarization map at visr —-10 km s~^, which is 
the same velocity map presented by Lai et al. (2003). In our case the visr —-10 km s~^ map 
does not present the same spatial distribution of polarization as the -8 km s~^ map, almost 
certainly due to the limited sensitivity of the polarization data. Table 1 gives a comparison 
between P.A. for both CO transitions. We can see that there is a consistent 90° difference 
between the transitions. 



2.5. Comparison between CO J = 1 — > and CO J = 2 — > 1 polarization 

Lai et al. (2003) presented a single velocity channel map of their line polarization 
observations (at visr — —10 km s~^); this is the velocity of the peak value in the Stokes I 
emission. Our CO J = 1 — > observations have a peak Stokes I at visr — —8 km s"^ (Figure 
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3 ), which also coincides with our most complete polarization spatial distribution. 

The Lai et al. (2003) map shows distributed CO emission that is more extensive than 
the region of the 1.3 mm continuum sources. We see a similar situation in our observations 
(Figure 3). We also see that our line polarization has better coverage over the CO emission 
than the Lai et al. (2003) observations (Figure 4). Our observations seem to indicate that 
with better signal to noise, it might be possible to achieve complete polarization coverage 
over the region of CO emission. The maximum polarized intensity coincides with the position 
of the continuum source, a fact which suggests the importance of a continuum radiation field 
in polarizing the line emission. This was studied in our numerical calculation. 

Comparing polarization vectors from Figure 3 and Figure 4 with the map presented by 
Lai et al. (2003), we found that the P. A. differ; in the central region in both maps most 
of the P.A. are orthogonal. This is a rather unexpected discovery. The Goldreich - Kylafis 
effect predicts that polarization can be either parallel or perpendicular to the magnetic field 
(Goldreich & Kylafis 1981a), but did not make a distinction between polarization of the 
same molecule at different transitions. This is because their calculation was done simulating 
a 2-level molecule. Deguchi & Watson (1984) extended the calculation to a multilevel 
molecule, but made no prediction about polarization direction in different transitions. Direct 
comparison of angles in Table 1 shows the 90° difference between both CO transitions. This 
is particularly true in the central part of the map. 

3. Polarization of molecular lines 

3.1. The Goldreich - Kylafis eff'ect 

The Goldreich - Kylafis effect is described in a series of papers (Goldreich & Kylafis 
1981a,b; Kylafis 1983a,b). They considered a molecule with only two (rotational) states 
having angular momenta 1 and 0, and investigated line formation and polarization in the 
presence of a magnetic field and anisotropic optical depths. Their original work predicted 
linear polarization up to about ten per cent. Motivated in part by an unsuccessful obser- 
vational survey for this polarization (Wannier et al. 1983), Deguchi & Watson (1984) 
extended the calculations to include a number of rotational states and radiative transitions. 
As a result, the predicted linear polarizations were reduced, typically by a factor of about 
two. We have further extended the multilevel calculations by including an external source 
term to represent emission from dust in a compact source. We have also used improved rates 
for collisional excitation, though this proved to be unimportant. We focus on calculations 
to understand the polarization angles of the CO J = 2 — > 1 and J = 1 — > transitions 
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observed in DR21(0H). 



3.2. Basic methods 

The formulation used here foUows Deguchi & Watson (1984). The radiative transfer 
equations are solved in the large velocity gradient or Sobolev approximation (hereafter LVG) , 
and in the regime where the line width Ai/, the Zeeman splitting g/ioB, and the natural line 
width 7 obey the relation Au » gHoB » 7. This strong inequality is easily satisfied here. 
The quantization z axis is along the magnetic field direction, and the quantum states are 
specified in the usual way by the total angular momentum J (there is no fine or hyperfine 
structure here) and by its projection M on the z axis. Under these conditions the radiative 
transfer equations for radiation associated with transitions between upper state J and lower 
state J' can be written as 



where || (_L) indicates the intensity of the radiation with a linear polarization parallel (per- 
pendicular) to the magnetic field. 

The difference with Goldreich & Kylafis (1981a) comes in expressing the opacity and 
source terms in a way to allow for arbitrary angular momenta. The opacity and source terms 
are thus written as hjmj'M' and Sjmj'M' such that Kjj, and Sjj, (where q stands for || or 

I^JJ' = \(t>{l^ - ^JJ') XI 'SJMJ'M' (3) 
AM=1 

c,± _^AM=l''^JMJ'M'SjMJ'M' /,^ 

-^JJ' - l^J 

Z^AM=1 I^JMJ'M' 



{ Sin^e ^ KJMJ'M' + ^COS^O ^ KJMJ'M' j (5) 
V AM=0 AM=1 ) 

II _ sin^O^^j^^QKjMJ'M'SjMJ'M' + ^COs'^O^^j^^'^KjMJ'M'SjMJ'M' 
Sm^O 2^AM=0 I^JMJ'M' + 2^0S^{U) 2^AM=1 I^JMJ'M' 
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where the symbol 9 represents the angle between the magnetic field and the direction of 
propagation and kjmj'm' and Sjmj'm' are defined in terms of Einstein A-coefficients and 
the populations tijm per magnetic substate 



2 

C ' 



K'JMJ'M' — TT- AjMJ'M'{nj'M' — ^jm) (7) 

87r ) 

^ _ hujj, njM , . 

iJJMJ'M' — 5 ■ \o) 

& nj'M' - njM 

Some equations here differ slightly from those in Deguchi & Watson (1984) because we are 
treating all magnetic substates explicitly. All statistical weight factors are one and are thus 
omitted for simplicity. 



Since the radiative transfer equations are functions of the populations, rate equations 
must be solved for these populations. We assume steady state , ^ = 0. The rate equations 
for the populations per magnetic substate are 

j<i j>i j 

where the Pij (which involve the radiation) are given in terms of the Einstein Ajmj'M' 
coefficients and the stimulated emission coefficients Rij. The indices i,j correspond to 
(J, M), (J', M'). The Cij are the coUisional excitation rates from state i to state j. Deguchi 
& Watson (1984) used the Cij given by Green & Chapman (1978). In our calculation 
we have updated the values of with the values provided by Flower (2001), which are 
given for a wide range of gas temperatures. We used a weighted average for the contribution 
of ortho and para hydrogen cross sections, specifically, Cij = 0.7 x Cf^^"" + 0.3 x Cfj^^". 
We considered temperatures up to Tgas = and angular momenta up to J = 9. The 

expressions for the Pij are 

Pij = Aij [ui + Rij{-nj + Ui)]. (1 0) 

The population equations constitute a nonlinear system, which we solve by numerical itera- 
tion. 

In the limit in which the macroscopic velocity differences of the gas in the cloud are 
greater than the thermal velocities of the molecules, the LVG approximation can be used to 
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simplify the solution for the radiative transfer equations. The LVG approximation expresses 
the intensity as a function of local variables. The intensity that emerges from the gas cloud 
and is detected by an observer is 

= {s% - f ) [1 - exp(-rj,,)] (11) 

where B represents the cosmic background radiation which is included in the calculation, 
and rjj, is the LVG optical depth for radiation with polarization q at the frequency of the 
Jj' transition, 

rjj, = 4,, X m (12) 
and where L{6) is the LVG characteristic scale length for which we adopt the form, 

L{e) = LQ/{asm^{e) + cos^(^)). (13) 

When a « 1,L{6) reflects a velocity gradient that is perpendicular to the magnetic 
field; for a >> 1, the velocity gradient is along the magnetic field (e.g., Deguchi & Watson 
(1984)). In the calculation, the constant Lq enters multiplied by the total number density of 
CO molecules. We vary this product to obtain the solutions as a function of optical depth 
that are presented in the Figures. 

The Rij in equation 10 involve the integral over frequency and over solid angle of the 
intensity within the cloud at the location of the CO, and are non-zero only for \ J — J'\ = 1, 
|M - M'l = and |M - M'l = 1 (Deguchi & Watson 1984), 

Q(-2 r HQ C II 

^2h7^J -^''""'^ J '^^'^^^ " (14) 



^Ahv^j^j ^'^'^^'^ " ^JJ'MA^) + <^osm]jm- (15) 
In the LVG approximation the integrals over frequency reduce to 



/ 



I%ct>{v - vjj.)dv = S%{1 - (3%) + ^^^^(5% (16) 
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in which the escape probabihty function is = [1 — Tjj,)]/Tjj,, and S(J1) represents 
the radiation that is incident from a compact, external continuum source. We express S(Q) 
as 



S{n) = (1 - e-^^)B,{Tsource) (17) 

for directions Q that arc subtended by the external source as viewed from the location of the 
CO. In other directions, S{Q) = 0. In equation 17, B,y{Tsource) is the Planck function, Tgource 
is the source temperature, and Tc is the continuum optical depth of the source at the relevant 
frequency. We adopted the parameterization of Mangum et al. (1991) for the combined 
emission from the dust core of DR21(0H), specifically, Tgource — 42 K, Tc=1.5[0.25/A[mm]]^ 
and an effective source angular size (radius) of 10" as seen from the earth. At the location of 
the CO, the angular size (radius) for the source is then approximately 1/2 radian (assuming 
0.3 pc for the separation between the CO and the source). Finally, the fractional polarization 
is expressed as 



— ^. (18) 

^± + ^11 



3.3. Calculation results 

The objective of this calculation is to understand why the the directions of the polar- 
izations for the CO J = 1 — and the CO J = 2 1 radiation arc perpendicular to 
one another in the DR21(0H) observations described in Section 2, with the J = 1 — > 
polarization parallel to that of the emission by the dust grains. We want to calculate the 
linearly polarized radiation emitted by CO for physical conditions that are representative 
for the DR21(0H) star forming region, and determine whether the observed directions and 
magnitudes of the polarization can be reproduced. The orthogonality of the polarizations is 
reflected by a sign difference in P in equation 18 . If P > 0, /_l dominates the emission and 
the polarization is perpendicular to the magnetic field. In the other case it is parallel. 

The physical conditions for our multilevel calculation are described below. These con- 
ditions are consistent with the assumption that the polarized emission comes from a cold 
envelope around the DR21(0H) continuum sources. We assumed that the physical conditions 
are the same in the gas for both CO transitions. 

To obtain a qualitative understanding of how the polarizations can be orthogonal for the 
two transitions, consider the J = 1 — > transition by itself, and first without any external 
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radiation. If the velocity gradients are smaller in directions parallel to the magnetic field 
than in directions perpendicular to the magnetic field, the optical depths for spectral lines 
will be the largest parallel to the field lines. The escape of radiation involved in de-exciting 
the upper J = 1 state will then be reduced more in directions along the field lines than in 
directions perpendicular to the field lines. This in turn leads to populations of the M = ±1 
substates that are larger than the populations of the M = substate because of the difference 
in the angular distributions of the de-exciting radiation. The angular distribution of the a 
radiation {\M — M] — 1 transitions) peaks in directions along the field fines whereas that of 
the TT radiation (|M — M] — transitions) peaks in directions perpendicular to the magnetic 
field. Hence, when the velocity gradients are smallest along the magnetic field, the rate for de- 
excitation of the M = ±1 substates is decreased more by trapping of the radiation than is the 
de-excitation rate of the M = substate . The excitation rate for all magnetic substates by 
collisions is the same. Under isotropic conditions, the populations of the magnetic substates 
arc equal. For any direction, the a and vr radiation that emerges from the radiative decays 
of these states will combine to give zero net polarization. However, when the M = ±1 states 
have larger populations, the a radiation will be relatively stronger in comparison with the 
TT radiation. Their contributions to the polarization will not then cancel and the combined 
radiation will be linearly polarized — in the direction of the polarization of the a radiation 
which is perpendicular to the magnetic field. 

Now consider by itself a compact, external source of radiation (i. e., no anisotropy in the 
optical depths as in the above paragraph) that is located in a direction perpendicular to the 
magnetic field. Because the angular distribution associated with tt transitions is larger for 
this direction than is the angular distribution of the a transitions, the J = 1, M = substate 
will be excited more rapidly by absorption of this radiation than will be the J = 1, M = ±1 
substates. The J = 1,M = substate would then become overpopulated in comparison 
with the other two magnetic substates and radiation that emerges from the gas would be 
polarized in the direction of tt radiation — parallel to the magnetic field. 

Though the reasoning is more complicated because of the additional substates and 
transitions, analogous conclusions apply for the polarization of the radiation associated with 
the J = 2 — > 1 and higher transitions of the CO molecule. 

When the above anisotropic velocity gradients and external source of radiation are 
both present, the direction of the polarization of the radiation that is emitted by the gas 
is determined by the relative importance of these two causes of anisotropy. Their relative 
importance is not, however, the same for the J = 2 — > 1 as for the J = 1 — > transition. The 
dust radiation from the core of DR21(0H) is much higher at the CO J = 2 — > 1 transition 
frequency than the CO J = 1 — > frequency, and thus the external radiation will play a 
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greater role for the polarization of the J = 2 — > 1 transition than it does for the J = 1 — > 
transition. The effect of the anisotropic velocity gradients is about the same for the two 
transitions since the collisional excitation rates are about the same for both transitions (the 
excitation energies for these states of CO arc less than kTgas)- The J = 2 — 1 radiation can 
then reflect the influence of the radiation and be polarized parallel to the magnetic field under 
conditions for which the J = 1 — > radiation reflects the influence of the velocity gradient 
and is polarized perpendicular to the magnetic field. The results presented below from the 
numerical calculations demonstrate that this orthogonality in the polarizations occurs for 
physical conditions that are likely for the gas that is emitting the polarized CO radiation 
from DR21(0H). Velocity gradients that are anisotropic and systematically larger in the 
direction perpendicular to the magnetic field occur in magnetohydrodynamic turbulence 
(see Watson et al. (2004)) 

Calculations were performed for gas temperatures between 15 and 50 K and for a range of 
densities. Dickel et al. (1978) obtained a Tgas—28 K; this temperature is used as a reference 
for our calculations. The velocity gradient must be largest in the direction perpendicular 
to the magnetic field according to the above reasoning. We thus focus on a = 0.1 in the 
expression for L{9), though we also present results for a=0.3 for comparison. 

Three curves for polarization versus LVG optical depth that correspond to calculations 
with densities = 25, 75 and 225 cm^'^ arc shown in Figure 5 for Tgas = 28 K. These 
densities are in agreement with what we might expect in a lower density envelope. Padin 
et al. (1989) obtained a density of = 5 x 10^ cm"'^, but this corresponds to the 
MMl continuum source and not to the DR21(0H) envelope. At the lowest optical depths, 
the external continuum source dominates in determining the direction of the polarization 
for both transitions whereas at the highest optical depths this radiation is largely excluded 
and the velocity gradient determines the direction of the polarization for both transitions. 
However, because the continuum emission from the core of DR21(0H) is much higher at the 
frequency of the J = 2 — > 1 transition than at the frequency of the J = 1 ^ 0, the external 
radiation remains more important to higher optical depths for the J = 2 ^ 1 transition 
than for the J = 1 ^ 0. Hence, a range of optical depths occurs over which the polarizations 
are orthogonal. This range of optical depths is mostly from log(r) ~ to —1 in Figure 
5. In this range of optical depths, the fractional polarizations typically are a few percent 
and are thus consistent with the observed magnitudes for the fractional polarizations of the 
two transitions. These small r values are consistent with our assumption that the polarized 
CO emission arises from a cold, low density gas, which in turn is compatible with the 
likely conditions for the DR21(0H) envelope. At the highest optical depths, the fractional 
polarization tends to zero in Figure 5 for both transitions as expected. The polarizations 
in the second panel in Figure 5 (for a— 0.3) provide an indication of the sensitivity of the 
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calculations to the magnitude of the velocity gradients. The polarizations for both transitions 
are presented only as a function of the CO J = 1 — > optical depth. However, the CO 
J — 2 ^ 1 optical depth is similar in our calculations. 

The curves in Figure 5 correspond to the densities that represent our best results. Our 
calculations show that for densities over nijj^SOO cm~^ the fractional polarization decreases 
dramatically for both lines, vanishing almost completely for densities over n/^^^lOOO cm~^. 
At the higher densities collisions are more rapid and the continuum source is no longer an 
important source of anisotropy. On the other hand, densities below nH2=20 cm^^ would 
require an envelope too big and far away from the continuum source, which would make the 
central continuum source ineffective as an additional source of anisotropy. 

In Figure 6, we demonstrate explicitly that the polarization is positive for both transi- 
tions at all optical depths in our calculations when there is no continuum source so that the 
velocity gradients determine the direction of the polarization. The fractional polarization 
also tends to zero at large and at small optical depths as expected, and as found by the 
previous investigators. 

In Figure 7, we present the results of our calculations for the other extreme; the limit in 
which the radiation from the continuum source dominates. We cannot eliminate collisions 
completely in our calculations. However, we can greatly reduce their influence by multiplying 
the collision rates by a small factor (10~^). As expected, the polarization is now negative for 
both lines (parallel to the magnetic field), is independent of density and converges to zero 
for higher optical depths. 

In our calculations, we also varied the gas temperature in order to explore the infiuence 
of uncertainty about the gas temperature in the envelope around DR21(0H). We found no 
significant variation in the polarizations for Tg^s between 15 and 50 K. What we found is 
that the range in r for which the polarizations are orthogonal shifts slightly in r; there is a 
displacement to lower r for higher temperatures. A molecular envelope might be heated by 
external UV photons or cosmic rays. This insensitivity to gas temperature shows a degree 
of robustness for the results of our calculations. We have also verified that the results are 
insensitive to the temperature of the dust by performing computations with Tgource up to 60 
K. 

The polarizations presented in Figures 5-7 are for the radiation that is emitted by the gas 
at an angle of 90° with respect to the magnetic field. We have verified that the polarization 
characteristics are similar for emission at angles from 90° to about 40°. At smaller angles, 
the polarized fraction decreases significantly — approaching zero at 0°, as expected. 
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3.4. Magnetic field strength and cloud support 



Our work has shown that the hne polarization will trace magnetic field geometry at 
lower densities than the dust polarization. However, comparing with the results of Lai et al. 
(2003), we believe that there is a correlation between the field traced by the dust (core) and 
the field traced by the line (envelope). Our dust polarization result at 3 mm shows similar 
P.A. to the 1.3 mm dust polarization of Lai et al. (2003). Also we now know why the hne 
polarizations are perpendicular at two different transitions, which makes the CO J = 1 — > 
polarization consistent with the CO J = 2 — > 1. At different frequencies the line polarization 
is either parallel or perpendicular to the dust polarization. The agreement in the geometries 
of the dust and line polarization suggests that the field in the core and the envelope are 
connected and not independent. 

Lai et al. (2003) estimated the magnetic field strength to be about 0.9 mC for MMl 
{nn, ~ 1 X 10^ cm.-^) and 1.3 mG for MM2 (nn^ ~ 2x 10^ cm"''^). Crutcher (1999a) measured 
a magnetic field in the line of sight Bios ~ -0.4 mG for MMl and Bios ~ -0.7 mG for MM2. 
Chandrasekhar & Fermi (1953) predicted a magnetic field dependence proportional to p^^^ 



The dispersions in the measured velocities and P.A. are similar for the J = 1 — > and 
J = 2 — > 1 lines. Therefore, we can extrapolate the values of the magnetic field measured in 
the core to the envelope following a p^^^ law. The molecular hydrogen number density for 
the cores is in the order of n^jj = 10^ cm~^, and taking a density of — 100 cm~^ for the 
envelope, we get a magnetic field in the envelope 100 times smaller than the core value of ~ 
1 mG (at ~ 1 x 10^ cm~^), or about 10 fiG. We can use this value to estimate the mass 
to magnetic fiux ratio; its critical value is given by (Mouschovias & Spitzer 1976) 



This equation can be expressed as a function of volume density in an spherical model, using 
the molecular hydrogen mass, equation 20, and expressing the magnetic field in G units, and 
the distance in cm. We arrived at 



BchF = Q^J^P 



A0 



(19) 
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in which R is the radius of the cloud, n is the number density of hydrogen molecules and B 
is the magnetic field. Using a density n =100 cm~^, B—IQ fxC and a radius of R=0.3 pc, 
we obtained a mass-to-flux ratio of 0.13 times the critical value, and thus highly subcritical. 
Hence, the envelope is supported against gravitational contraction by the magnetic field. The 
value of R represents an angular distance of 0.3 pc of distance for DR21(0H). This radius 
may be underestimated due to resolution of more extended structure by the interferometer. 
However, single dish maps (Wilson & Mauersberger 1990) suggest that R < 1 pc. 

4. Summctry and conclusions 

DR21(0H) was observed in 3 mm dust continuum and CO J = 1 — > line emission. 
Comparing our line observations with Lai et al. (2003), we observed a consistent difference 
in polarization direction between the CO J = 2 ^ 1 and CO J = 1 ^ lines; they are 
perpendicular to each other over the central region (which coincides with the position of 
the MMl and MM2 continuum sources). We developed a code based on the Deguchi & 
Watson (1984) calculation in order to solve the radiative transfer equations and calculate 
the fractional hnear polarization for different transitions of the CO molecule. We found that 
the presence of a small continuum source will likely produce an increase in the anisotropy of 
the radiation field over the CO gas. An anisotropic radiation field will unevenly populate the 
magnetic sub-levels of the CO molecule, which will produce linearly polarized emission. We 
showed that the presence of a compact continuum source will produce a gradual change in the 
fractional polarization, with a sign change, as a function of optical depth. This would explain 
the orthogonal orientations of the CO polarization in different transitions. In the particular 
case of DR21(0H), the physical conditions that are consistent with the polarization data 
correspond to a hot continuum source {Tdust ^ 42 K) and a CO gas temperature of Tgas ~ 
28 K, for a uh^ ~100 cm-^. 

Line polarization observations can become a powerful tool to constrain the magnetic field 
geometry. We have seen that a clear understanding of the physical conditions in a molecular 
cloud will give more accurate information about the fractional polarization sign. A good 
knowledge of the sources of anisotropy in a cloud will help to understand and constraint the 
geometry in line polarization observations at multiple frequencies. 

We also found that line polarization traces low density material; this is important in 
order to connect the magnetic field geometries at different densities. In the particular case of 
DR21(0H) we believe that there may be a correlation between the field traced by the dust 



-16- 



at densities of 10^ cm~^ and the field traced by line polarization at 10^ - 10^ cm~^. However, 
we were only able to compare our calculation with polarization observations of CO molecular 
transitions. We believe that polarization information from additional molecules may probe 
the magnetic field geometry at intermediate densities giving a more accurate picture of the 
field morphology in star forming regions. 

Having information about the magnetic field at different densities is useful to test 
the gravitational state of equilibrium for the cloud. Polarization mapping can allow the 
Chandrasekhar- Fermi method to be used to obtain an estimate of the magnetic field strength 
and test whether the envelope is magnetically supported or not. The envelope of DR21(0H) 
appears to be magnetically supported. 

This research has been supported partially by NSF grants AST 02-05810 and AST 
99-88104. 
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Table 1. Position angle for CO J = 2 — > 1 and CO J = 1 — > 1 at velocity channel map 
V — -10 km s~^. Data was interpolated at a tolerance of 0.2" that corresponds to 
approximate 2.7 x 10~^ pc using a distance to DR21(0H) of 3 kpc. 



Offsets in arcsec 




Offsets in arcsec 




(1.0,2.0) 


78±5 


(1.0,2.0) 


3.6±9.2 


(-0.6,2.0) 


77±5 


(-0.6,2.0) 


4±9.6 


(4.2,4.0) 


87±7 


(4.0,4.0) 


-0.2±8.9 


(2.6,4.0) 


85±6 


(2.5,4.0) 


-5.2±6.2 


(1.0,4.0) 


83±7 


(1.0,4.0) 


-6.1±5.3 


(-0.6,4.0) 


84±7 


(-0.5,4.0) 


-4.6±6.2 


(9.0,6.0) 


-60±5 


(9.0,6.0) 


15.3±9.7 


(7.4,6.0) 


-70±6 


(7.5,6.0) 


13.0±7.6 


(5.8,6.0) 


-85±8 


(6.0,6.0) 


6.7±6.9 


(4.2,6.0) 


88±8 


(4.2,6.0) 


-3.5±5.5 


(2.6,6.0) 


87±8 


(2.5,6.0) 


-9.3±4.1 


(9.0,8.0) 


-50±7 


(9.0,8.0) 


18.4±9.6 


(7.4,8.0) 


-62±8 


(7.5,8.0) 


14.1±7.6 


(5.8,8.0) 


-76±9 


(6.0,8.0) 


8.0±6.6 
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Fig. 1. — Polarization map of DR21(0H) at 3 mm. The contours represent Stokes I at -0.02, 
0.02, 0.04, 0.07, 0.09, 0.11, 0.14, 0.16 and 0.18 Jy beam-^ The pixel gray scale show 3(t 
polarized intensity (y^Q^ + U^) for the dust continuum emission also in Jy beam~^, while the 
black vectors are the dust fractional polarization. The length of the vectors are proportional 
to fractional polarization with the length of the bar at the top-right corresponding to a 
fractional polarization of 0.33. 
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Fig. 2. — The spatially averaged Stokes I spectrum from the CO J = 1 — > emission. 
We observe a peak of 5.6 Jy beam~^ at visr — —8 km s~^, while the minimum is located 

at visr — —2 km s~^. {upper panel). Plotted over is the averaged P.A. over the Stokes I 
spectrum. Most of the velocity channels have P.A. that are orthogonal to the Lai et al. 

rrtnno\ 1-1- /; A T-ii „ • j r ..• i l,, ^-.^-i-; ; 1„-i-i-„J 0-i-„1 
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RA {J2000) 

Fig. 3. — The CO J = 1 — > polarization emission from DR21(0H) at visr km s~^. The 
contours represent Stokes I at -6.15, -4.61, -3.08, -1.54, 1.54, 3.08, 4.61, 6.15, 7.7, 9.23, 10.77 
and 12.31 Jy beam" . The gray scale represents 3a polarized emission (v/Q2 + [/2) also m 
Jy beam^^. The vectors are a representation of the fractional polarization, the length of the 
bar at the top-right corresponds to a fractional polarization of 0.086. The most widespread 
polarization pattern occurs at this velocity, which coincides with the peak of the CO emission 
(Stokes I image). 
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Fig. 4. — The CO J = 1 ^ polarization emission from DR21(0H) at vis,- =-10 km s~^. 
The contours represent Stokes I at -2.73, -1.82, -0.91, 0.91, 1.82, 2.73, 3.64, 4.56, 5.47, 6.38 
and 7.3 Jy beam~^. The gray scale represent 3(T polarized emission 

WQ'^ + U^) also in Jy 

bearn"^. The vectors are a representation of the fractional polarization, the length of the 
bar at the top-right corresponds to a fractional polarization of 0.064. This map represent 
the same velocity channel presented by Lai et al. (2003). Comparing with Lai et al. (2003) 
at the same velocity channel we see that most of the polarization is perpendicular to the CO 
J = 2 — > 1 polarization. 
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Fig. 5. — Fractional polarization P calculated for the CO transitions J = 1 — > and 
J = 2 — > 1 as a function of the LVG optical depth for the CO J = 1 — > transition. Here, 



gas 



=28 K and T„ 



=42 K. Each pair of curves represents calculations for a different 



density of molecular hydrogen — 25 cm , 75 cm and 225 cm . In the upper panel, the 
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Fig. 6. — Fractional polarizations P calculated for the CO transitions J = 1 ^ and 
J = 2 — > 1 are shown as a function of r for the molecular hydrogen densities 25, 75 and 
225 cm~^ when there is no compact, external source of continuum radiation. Here, the 
direction of the polarizations is determined by the direction of the velocity gradient and is 
perpendicular to the direction of the magnetic field. Note that the polarizations of the CO 
J = 1 — > and J = 2 ^ 1 transitions are never orthogonal (i. e. they always have the same 
sign.) 



-26- 




-0.4 I 1 1 ' ' ' 1 

-1 -0.5 0.5 1 1.5 2 

log('c) 

Fig. 7. — Fractional polarization P for the CO transitions J = 1 ^ and J = 2 — 1 are 
shown as a function of r when the collision rates used in Figure 5 are reduced by multiplying 
by a factor 10^^. For the relevant hydrogen densities (25 - 225 cm~^), the external continuum 
source dominates in determining the direction of the polarization for both transitions and the 
polarizations do not depend on density. For both spectral hnes, the polarization is parallel 
to the magnetic field (negative sign) as expected, and converges to zero for large values of r. 



